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Introduction 

The  realization  of  efficient,  reliable  picosecond 
closing  switches  will  make  possible  high  gradient  (1- 
3  GeV)  linacs  with  pulsed  accelerating  electrode 
structures.  Recently  a  promising  candidate  for 
picosecond,  high  voltage  switching,  the  gas  avalanche 
switch,  has  been  proposed. !>2  The  medium  in  this 
switch  is  high  pressure  (10-800  atm)  gas.  An  avalanche 
discharge  is  initiated  between  pulse-charged  high 
voltage  electrodes  by  multiphoton  ionization  from  a 
picosecond  order  laser  pulse.  The  laser-initiated 
electrons  avalanche  toward  the  anode,  causing  the 
applied  voltage  to  collapse  in  picoseconds,  long  before 
the  hot  current  channel  formation  of  a  conventional 
spark  gap. 


permittivity,  n  is  the  electron  density,  e  is  -1.6  x 
10-19  coui j  vx  is  the  electron  velocity  in  the  x 
direction,  Vy  is  the  electron  velocity  in  the  y 
direction  and  pg  is  the  permeability. 

The  electron  fluid  equations  for  continuity, 
momentum,  and  energy  conservation  are 
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Several  versions  of  the  gas  avalanche  switch  may 
be  conceived.  A  parallel  plate  capacitor  version 
consists  of  a  high  pressure  gas  confined  between  two 
parallel  metal  plates.  An  analysis  of  the  formation  or 
voltage  delay  and  voltage  collapse  times  for  this 
geometry  has  been  done  by  Villa  and  Cassell1*2  by 
numerical  integration  of  a  zero-dimensional  (QD) 
electric  circuit  equation,  which  includes  a  resistively 
limited  charge  feed.  This  analysis  predicts  voltage 
collapse  times  of  the  order  of  several  picoseconds  for 
N2  and  air  and  1  psec  or  less  for  Ar. 

The  Two-Di mens ional  E le ctromagnetic 
Electron  Fluid  Avalanche  Code 

For  confirmation  of  the  encouraging  results  of  the 
circuit  code,  we  have  written  a  two-dimensional  (ZD) 
code,  which  simultaneously  and  sel fconsistently  solves 
Maxwell's  equations  for  transverse  magnetic  (TM)  field 
modes  in  a  gas  between  two  parallel  perfect  conductors 
and  a  set  of  electron  fluid  conservation  equations.  A 
rectilinear  perfect  conductor  may  be  placed  between  the 
pi  ates . 

For  a  free  space  dielectric  between  the  plates, 
the  relevant  Maxwell's  equations  in  the  MKS  system  are 
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where  Ex  is  the  electric  field  in  the  x  direction,  Ey 
is  the  electric  field  in  the  y  direction,  H,  is  the  1 
magnetic  field  in  the  z  direction  (the  unmodeled 
direction  perpendi cul ar  to  x  and  y),  eg  is  the 


^  =  ne(ExVx+  EyVy)  -  nVu(U-U0)  -  "Vi  (7) 

where  v.j  is  the  ionization  collision  frequency,  m  is 
the  electron  mass,  vm  is  electron-neutral  molecule 
momentum  transfer  frequency,  U  is  the  average  electron 
kinetic  energy,  vu  is  electron-neutral  energy  transfer 
frequency,  Up  is  the  average  neutral  energy,  and  is 
the  neutral  ionization  potential.  The  equations  are 
specialized  to  air  by  taking  Uo  -  0.025  eV  and 

=  14  eV.  It  must  be  noted  that  convection,  pressure 
gradients,  heat  flow,  and  magnetic  forces  are 
neglected.  Electron  losses  due  to  atomic  and  molecular 
processes,  such  as  attachment  and  recombination ,  are 
ignored.  Electron  generation  from  photoionization  is 
also  ignored. 

The  collision  frequencies  for  air  are  those  used 
previously  in  microwave  breakdown  calculations.2  These 
frequencies  are  assumed  to  scale  linearly  with  neutral 
pressure.  In  addition  electrode  pulse  charging  is 
assumed  to  be  instantaneous  so  that  the  very  real 
experimental  problem  of  premature  breakdown  due  to 
effects,  such  as  field  emission,  is  ignored.  The 
initial  laser  deposition  of  electrons  is  also  assumed 
instantaneous.  Thus  the  details  of  modeling  the  laser 
ionization  are  presently  conveniently  bypassed. 

The  space  between  the  horizontal  parallel  plates 
is  partitioned  with  a  rectangular  mesh.  Solution 
points  for  the  electromagnetic  and  fluid  variables  are 
sited  to  obtain  a  staggered  grid,  similar  to  those  in 
explicit,  2D,  finite  difference,  electromagnetic 
scattering  codes.  Ex  is  solved  at  even  points  along 
odd  horizontal  (xz)  planes.  The  lowest  of  these  planes 
coincides  with  the  bottom  metal  plate;  whereas,  the 
uppermost  plane  coincides  with  the  top  capacitor 
plate.  Ey  is  solved  at  odd  points  along  even 
horizontal  planes;  whereas,  Hz  is  solved  at  even  points 
on  the  same  planes.  The  four  fluid  variables  are 


* Th i s  work  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the  Lawrence  Livermore  National 
Laboratory  under  Contract  W-7405-Eng-48 . 
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solved  at  points  where  Hz  is  defined.  The  spatial 
derivatives  in  Maxwell's  equations  are  replaced  with 
central  finite  differences.  Seven  ordinary 
differential  equations  in  time  result  for  each  grid 
cell  of  four  grid  points.  The  global  set  of  equations 
for  all  the  grid  cells  is  solved  by  time  integration 
with  the  block-iterative,  stiff,  implicit,  variable 
step  solver  GEARBI.4 

Modeling  Parallel  Plate  Capacitor  Switches 

The  2D  avalanche  code  was  first  applied  to  a 
parallel  plate  capacitor  with  a  spacing  or  gap  of 
100  p  and  an  extent  in  the  x  direction  of  247.5  u.  An 
initial  electron  number  of  10^  (n  =  4  x  10^1  m_3)  Was 
uniformly  distributed  for  an  assumed  depth  of  1  m  in 
the  z  direction.  A  voltage  of  266  kV  was  applied  to 
the  top  plate  by  initialization  of  a  spatially  uniform 
Ey.  Boundary  conditions  were  applied  at  x  =  0  and 
247.5  p,  which  approximated  the  infinite  extension  of 
the  parallel  plate  capacitor.  An  air  pressure  of 
350  atm  was  specified,  and  the  simulation  was  started. 

The  resulting  plate  voltage  waveform  is  shown  in 
Figure  1,  labeled  "2D  Code."  The  voltage  across  the 
plates  stays  almost  constant  for  about  3  psec,  then 
decays  exponentially.  The  voltage  collapse  time  is 
about  2  psec.  The  result  of  a  circuit  code  calculation 
by  Cassell  and  V i 1 1  a 2  is  shown  as  the  curve  labeled 
"Circuit  Code."  For  this  calculation,  the  first 
Townsend  coefficient  a  was  set  to  2/3  of  its  reported 
experimental  value. 


Figure  1.  Comparison  of  voltage  waveforms  at  266  kV  and 
350  atm.  Alpha  is  the  first  Townsend  coefficient  used  in 
the  circuit  code  calculations  and  is  equal  to  0.67  times  the 
reported  value. 


This  adjustment  is  supported  by  Figure  2,  which 
shows  the  ratio  of  the  pressure-normal i zed  first 
Townsend  coefficient  a/p  from  the  2D  Code  to  that  from 
experimental  data. 5  The  quantity  p  is  the  neutral  gas 
pressure  in  Torr.  The  stored  electric  field  energy  is 
also  shown  as  a  percentage  of  the  initial  electrostatic 
field  energy.  Over  most  of  the  decay  of  the  electric 
field  energy,  a  from  the  2D  code  is  about  0.6-0. 7  of 
the  experimental  value. 


E/p  (V/cm-Torr) 

Figure  2.  Ratio  of  p-  for  the  code  to  measured  p-. 

Figure  3  shows  the  average  electron  velocity  from 
the  2D  code  versus  E/p  for  the  range  0.1  -  16.5  V/(cm- 
Torr).  The  electron  drift  velocity  measured  by 
Bradbury  and  Nielson^  is  also  shown.  The  agreement  of 
the  two  curves  is  excellent. 


Figure  3.  Comparison  of  electron  drift  velocities. 


Figure  4  shows  the  change  of  the  capacitor  voltage 
waveform  with  pressure  at  constant  initial  voltage  and 
electron  number.  The  initial  E/p  is  shown  with  each 
pressure.  As  the  air  pressure  increases  and  the 
initial  E/p  consequently  decreases,  both  the  delay  time 
and  the  voltage  collapse  time  increase. 

Further  calculations  show  picosecond  scale  delay 
and  collapse  times  on  the  same  100  p  plane  gap  at  lower 
and,  hence  more  experimentally  appealing,  charging 
voltages.  At  106  initial  electrons/m  of  plate  depth, 

80  kV,  and  27.2  atm,  the  voltage  delay  time  is 

2.2  psec;  the  decay  time  is  1.8  psec.  At  the  same 

electron  number,  62  kV,  and  17.0  atm,  the  delay  time  is 

2.3  psec;  the  decay  time  is  1.7  psec.  At  43  kV  and 

10.2  atm,  the  delay  time  is  2.8  psec,  and  the  decay 

time  is  2.1  psec.  These  conditions  of  voltage  and 
pressure  are  relevant  because,  in  recent  preliminary 
experimental  work  with  M2,  a  pulse-charged  102  p  curved 
to  xnife  edge  gap  has  withheld  these  voltages  for 

6  nsec  with  an  irregular  ac  pulse. 
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Time  (psec) 

Figure  4.  Voltage  waveforms  versus  air  pressure. 


Mode  1 i ng  TEM  Wave  Generat i ng  Accelerator  Switches 

A  more  sophisticated  gas  avalanche  switch, 
suitable  for  a  pulsed  power  linac,  can  be  made  by 
placing  a  properly  dimensioned  center  electrode  between 
the  parallel  plates  and  charging  it  to  a  suitable  high 
voltage.  A  simple  model  of  such  a  linac  switch  is 
shown  in  Figure  5.  The  parallel  plates  are  initially 
grounded;  the  air  between  them  is  at  700  atm.  The  gap 
between  the  center  electrode  and  the  bottom  plate  is 
uniformly  filled  with  1.7  x  1013  electrons  for  an 
assumed  depth  of  1  m  into  the  figure.  For  a  center 
electrode  voltage  of  3.889  MV,  this  configuration 
produces  electromagnetic  waves,  which  are  predominantly 


TEM  modes  and  propagate  outward  from  the  center 
electrode  to  the  left  and  right  boundaries. 

The  initial  electric  field  components  Ex  and  Ey 
are  obtained  from  the  finite  element,  electrostatic 
code  STAT2D  and  used  as  initial  conditions  for  the 
electron  avalanche  code.  After  the  simulation  starts, 
the  voltage  between  the  center  electrode  and  the  lower 
plate  collapses  within  1  psec.  Avalanche-generated  Ey 
waves  at  10  psec  are  shown  in  Figure  6.  The  wave 


Time=1 0.000  1 0'1!  sec 


crests  in  the  near  Ey-x  plane,  which  corresponds  to  the 
lower  parallel  plate,  will  hit  the  x  boundaries  at 
about  12  psec.  The  wave  crests  in  the  far  Ey-x  plane, 
which  corresponds  to  the  top  plate,  lag  behind  those  in 
the  near  Ey-x  plane.  The  bipolar  spikes  below  the 
center  electrode  are  from  Hz  noise  introduced  from 
inaccuracies  in  the  extrapolated  STAT2D  electric  fields 
at  the  center  electrode  surface.  In  addition  to  the 
strong  Ey  waves,  weak  Ex  waves  propagate  from  the 
center  electrode. 

The  voltage  waveform  of  the  top  plate  with  respect 
to  the  lower  plate  at  the  right  boundary  is  shown  as 
the  solid  curve  in  Figure  7,  where  the  rise  time  is 
less  than  2  psec.  The  first  peak  is  from  the  crest  at 
the  lower  plate;  the  second  peak  is  from  the  lagging 
crest  at  the  top  plate.  The  voltage  full  width  at  half 


Figure  7.  Plate  voltage  waveform  at  the  right  boundary. 
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maximum  (fwhrn)  is  about  8  psec.  The  peak  voltage  of 
4.769  MV  is  a  22.6  %  enhancement  of  the  initial  center 
electrode  voltage.  This  enhancement  is  most  likely  due 
to  the  transmission  properties  of  the  transmission  line 
transitions  at  the  sides  of  the  center  electrode,  where 
two  short,  separated  narrow  lines  join  a  wider  line. 

At  29.5  psec,  422  J  or  71.4  %  of  the  initial 
electrostatic  energy  has  been  transported  across  the 
boundaries . 

When  1.2  x  10^  initial  electrons  are  evenly 
concentrated  into  a  50  p  wide  column  centered  under  the 
center  electrode,  a  peak  right  boundary  voltage  of 
5.378  MV  and  the  dashed  waveform  in  Figure  7  result. 
This  is  a  38.3  %  voltage  enhancement.  In  this  case 
490  J  or  82.9  %  of  the  initial  energy  has  been 
transported  out  of  the  computational  space  by 

29.5  psec. 

Further  simulations  predict  that,  without 
convection  or  photoionization,  the  center  electrode  to 
lower  plate  space  must  be  completely  bridged  by  initial 
electrons  for  maximal  voltage  pulse  generation  and 
energy  transport.  If  1.7  x  i0i3  initial  electrons  are 
concentrated  uniformly  in  the  two  thirds  of  that  space 
closest  to  the  center  electrode,  then  the  peak  pulse 
voltage  is  3.15  MV,  and  the  energy  transport  is  20  %. 

If  the  same  initial  electrons  are  concentrated 
uniformly  in  the  one  third  of  that  space  closest  to  the 
center  electrode,  then  the  peak  pulse  voltage  is  2.11 
MV  and  the  energy  transport  is  9.4  %. 

The  previously  mentioned  experimental  work  has 
shown  encouraging  voltage  holdoff  in  N2  at  80  kV  and 
27.2  atm.  These  conditions  represent  an  initial  plane 
gap  E/p  of  379  V/(cm-Torr).  At  300  k V  and  27.2  atm, 
the  initial  nominal  center  electrode  E/p  on  the 
accelerator  switch  in  Figure  5  is  386  V/(cm-Torr)  with 
a  600  p  gap.  A  simulation  at  these  conditions,  where 
there  is  greater  expectation  for  sufficient  voltage 
holdoff  time,  produces  a  boundary  voltage  pulse  of  2.5 
psec  risetime  and  about  9  psec  fwhrn.  The  peak  pulse 
voltage  is  301  kV  and  the  energy  transport  is  55.4  %  at 

36.5  psec. 
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Conclusions 

The  parallel  plate,  2D  electron  fluid  avalanche 
code  predicts  voltage  collapse  times  of  1-10  psec  for 
air  pressures  of  80-500  atm.  The  collapse  times  agree 
well  with  0D  circuit  code  results.  The  electron  drift 
velocity  from  the  2D  code  agrees  well  with  experimental 
results  for  0.1  s  E/p  <  16  V/(cm-Torr).  The  a/p  from 
the  2D  code  is  in  good  agreement  with  experimental 
results  for  46  ^  E/p  S  100  V/(cm-Torr).  Simulations  at 
lower  voltages  and  pressures,  for  which  there  are  some 
promising  preliminary  experimental  results  for  voltage 
holdoff  at  the  same  nominal  electrode  gap,  also  give 
picosecond  scale  (1.7-2. 1)  voltage  collapse  times. 

A  linac  switch  structure  with  a  parallel  plate 
spacing  of  2  mm,  a  pressure  of  700  atm  air,  an  initial 
nominal  center  electrode  E/p  of  200  V/(cm-Torr)  and  a 
charging  voltage  of  3.889  MV  produces  voltage  collapse 
in  less  than  1  psec  and  TEM  wave  generation. 

Transmitted  voltage  pulse  risetimes  of  2  psec,  voltage 
pulse  enhancements  of  23-38  %,  and  energy  transports  of 
71-83  %  of  the  initial  electric  field  energy  are 
predicted.  Lower  voltages  (300  kV)  and  pressures 
(27  atm)  on  this  very  same  switch  structure  give 
similar  results . 
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